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O ■ Abstract 



We study junctions consisting of confining strings in M = 1 su- 
C^i! persymmetric large N gauge theories by means of the gauge/gravity 

correspondence. We realize these junctions as D-brane configura- 
tions in infrared geometries of the Klebanov-Strassler (KS) and the 
Maldacena-Niinez (MN) solutions. After discussing kinematics asso- 
ciated with the balance of tensions, we compute the energies of baryon 
vertices numerically. In the KS background, baryon vertices give nega- 
tive contributions to the energies. The results for the MN background 
strongly suggest that the energies of baryon vertices exactly vanish, 
as in the case of supersymmetric (p, g)-string junctions. We find that 
brane configurations in the MN background have a property similar to 
the holomorphy of the M-theory realization of (p, g)-string junctions. 
With the help of this property, we analytically prove the vanishing of 
the energies of baryon vertices in the MN background. 



* E-mail: imamuraOhep-th . phys . s . u-tokyo .ac.jp 



1 Introduction 

The gauge/gravity correspondence^ El 0] relates five- dimensional gravitational 
backgrounds to four- dimensional field theories on the boundaries of the five- 
dimensional spacetimes. The extra dimension r is related to the energy scales 
of the field theories through red-shift (warp) factors depending on r. The five- 
dimensional spacetimes are in general accompanied by compact internal spaces. 
The structures of these higher- dimensional spacetimes reflect non-perturbative 
properties of their dual field theories. Many gravity duals of various field theories 
have been constructed in the context of string theory as near horizon geometries 
of branes on which gauge theories are realized. 

Various kinds of particles in boundary field theories are identified with a 
number of objects in string theory. These particles have been investigated with 
respect to the above-mentioned duality. Specifically, the spectra of glueballsjUEl 
El III , mesons|BJ El EH EJJ EE2 ECU El IIS! and (di-)baryons[IEl 03 03 EH1 ED] and 
interactions among themj^DESE] have been studied with a variety of methods. 
Bound states of massive adjoint or bifundamental particles are investigated in 
Refs. [HI El E3 EI3 EH EH Ei- In Ref. mi, the duality is used to account for 
the extremely narrow decay width of pentaquark barvons pJH 132 j. 

Hadrons in SU(iV) gauge theories are bound states of (anti-)quarks belonging 
to the ( ant i-) fundamental representation. In order to introduce such quarks and 
antiquarks, we need to place flavor D-branes in the dual gravity background |33|. 
Quarks are identified with endpoints of open strings on the flavor branes, and 
hadrons are constructed by connecting them. For example, a meson consisting 
of a quark and an antiquark is realized as an open string stretched between two 
flavor branes on the gravity side, as depicted in Fig. E (a). 



Figure 1: (a) A meson configuration and (b) a baryon configuration. 

Among the many states in a meson spectrum, some low-lying states of pseudo- 
scalar and vector mesons are identified with string modes representing fluctua- 
tions of the flavor branes and gauge fields on them. The masses of such modes 
are examined in Refs. El] , and the existence of localized modes with a mass 
gap is found in those works. 

By contrast, Regge trajectories of higher-spin mesons can be treated as semi- 
classical spinning strings in curved backgrounds |Ul ITS] . In general, it is difficult to 
determine the motion of strings completely, and the Born-Oppenheimer approxi- 
mation is often used to simplify such problems. In this method, a quark-antiquark 
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(b) 
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potential is computed as the energy of a string with fixed endpointsj^H ESI, an d 
then, the motion of the endpoints (the quark and the antiquark) is obtained by 
solving the resulting potential problem. If the distance L between a quark and 
an antiquark is large, the L dependence of the potential energy can in general be 
expanded as 

E q q = TL + M q + Mg + 0{L~ l ). (1) 

The coefficient T of the first term here is the tension of the confining string be- 
tween the quark and the antiquark. From the gravity point of view, the first 
term is interpreted as the contribution from the tension of a fundamental string 
at r = 0, because fundamental strings tend to approach r = as a result of the 
gravitational force. Thus, T is obtained as the product of the redshift (warp) fac- 
tor and the proper tension of fundamental strings. In general, confining strings 
can be bound states of elementary confining strings. The tension of such strings 
depends non-linearly on the number of constituent elementary strings. This prop- 
erty is reproduced in the gravity description [SHI by taking account of the effect 
of the transition of fundamental strings to D-branes expanded in internal spaces 
by Myers' effect p?Tj . The second and third terms in (0), which are independent 
of L, are interpreted as the self-energies of the quark and antiquark. These can 
be evaluated by analyzing the catenary profile of a string near flavor D-branes. 

A baryon is constructed by connecting flavor D-branes with N open funda- 
mental strings of the same orientation. In order to join such strings without 
violating the conservation law of the fundamental string charge, we have to in- 
troduce a baryon vertex at the joint ^H] E] (see Fig. G] (b)). This baryon vertex 
is a certain D-brane wrapped around a non-trivial cycle in the internal space. In 
general, strings stretched between a baryon vertex and a flavor D-brane could be 
bound, and the number of confining strings meeting at a vertex can be smaller 
than N. 

Let us consider the potential among quarks in a baryon in the spirit of the 
Born-Oppenheimer approximation. [T7j If the length of each confining string is 
large, the energy of the baryon configuration is given by 

-^baryon = TjLj + ^ M, + -^vertex + ^(L' 1 ), (2) 

i i 

where Li and T\ are the length and the tension of the i-th confining string, and 
Mi represents the self-energy of the quark at the end of the i-th confining string. 
The third term, E vcvtcx , represents a new contribution which is absent in meson 
energies. Naively, this contribution can be evaluated by multiplying the tension 
T of the D-brane composing the baryon vertex, the volume V of the cycle around 
which the D-brane is wrapped, and the warp factor. Although this does not give 
an exact result, because of the deformation of the brane caused by the string 
tension [TBI EH], we expect that a vertex gives some positive contribution of the 
order of TV x (warp factor). The purpose of this work is to compute this quantity 
using both numerical and analytical methods. 
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This paper is organized as follows. In the next section, we describe the context 
we consider and clearly define the problem to be solved. The junctions are realized 
as D3-brane configurations embedded in a six-dimensional geometry. In £0 we 
briefly review the result of Ref. [SB] concerning the confining string tension and 
study the kinematics of junctions associated with the tension balance. In §3J as 
preparation for the following computations, we rewrite brane configurations as 
two-dimensional surfaces in a certain four-dimensional target space by assuming 
a certain symmetry. The results of numerical computations are reported in £0 
and we present analytic results in §B1 We conclude in §3 I n the appendix, we 
briefly explain how the electric fields on D3-branes are treated in the numerical 
computations. 

2 D3-branes in the infrared geometry 

We use two different classical solutions in type IIB supergravity as background 
spacetimes dual to M = 1 supersymmetric confining gauge theories. One is the 
Maldacena-Niinez (MN) solution (SHUSH], an d the other is the Klebanov-Strassler 
(KS) solution 0D| El E2|- 

The MN solution is the near horizon geometry of N coincident D5-branes 
wrapped around a two-cycle in a non-compact Calabi-Yau manifold. This is 
believed to be the gravity dual to the Af = 1 SU(iV) pure Yang-Mills theory. 
This theory has only two parameters, the size of the gauge group, N, and the 
confinement scale, Aqcd- 

The KS solution can be interpreted as the near horizon geometry of N + k 
coincident D5-branes and k anti-D5-branes wrapped around a two-cycle. The 
corresponding boundary field theory is an M = 1 SU(iV+ k) x SU(/c) gauge theory 
with bifundamental chiral multiplets. The number k depends on the energy scale. 
As the energy decreases, k also decreases due to cascade phenomenon |41j . In the 
brane picture, this can be regarded as a kind of pair annihilation of D5-branes 
and anti-D5-branes. At very low energy, all the anti-D5-branes disappear with 
the same number of D5-branes, leaving N D5-branes. Thus, the boundary theory 
of the KS solution behaves like the M = 1 SU(A r ) supersymmetric Yang-Mills 
theory in the low energy limit. In addition to N and Aqcd ; this theory includes 
the parameter g, defined by g~ 2 = g x 2 + g% 2 , where g\ and g 2 are the gauge 
coupling constants of SU(A r + k) and SU(fc), respectively. Although g x and g 2 
run with the energy scale individually, g is a renormalization invariant parameter. 
This parameter is related to the string coupling constant g str , which is constant 
in the KS solution, as 

g 2 = 4irg str . (3) 

Corresponding to the similarity of the types of low-energy behavior of the 
two gauge theories, the two classical solutions have similar structure. They are 
warped products of a five- dimensional manifold with coordinates (x M ,r) and an 
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internal space with the topology S 2 x S 3 . S 2 shrinks to a point at r = 0, while 
S 3 possesses non-zero size everywhere. The infrared (IR) dynamics of the field 
theories are reflected by the structure near the centers of the classical solutions. 
For both the MN and the KS solutions, the metric of the r = subspace relevant 
to the IR dynamics is given by 

ds 2 = R^^dy^dy" + dQl), (fji,u = 0,1, 2, 3) (4) 

where dVl\ is the metric of the unit 3-sphere. We refer to this as IR geometry. 
The size N of the low-energy gauge group SU(iV) of the boundary gauge theories 
is determined by the R-R 3- form flux flowing through S 3 . In other words, the 
integral of the R-R 3-form field strength G 3 over the non-trivial 3 cycle gives TV: 



N = G 3 . (5) 

The boundary theories of the MN and KS solutions include only adjoint and 
bifundamental fields. In order to introduce fundamental quarks into these theo- 
ries, we need flavor D-branes. We can use D5jl21II2], D7 [321 CHI, or D9-branes|4^] 
wrapped around appropriate cycles as flavor branes. Hadrons are constructed by 
connecting these branes with fundamental strings joined by baryon vertices. We 
assume that the lengths of the strings are large and focus on the vicinity of baryon 
vertices, which are located at the bottom (r = 0) of the classical solutions, due 
to the gravitational force. We do not discuss the effect of the endpoints of strings 
on the flavor branes, and the arguments given in this paper are independent of 
the choice of the flavor branes. 

In the cases of the MN and KS backgrounds, baryon vertices are D3-branes 
wrapped around S 3 EHl IS] ; and fundamental strings in the IR geometry are 
expanded to D3-brane tubes with section S 2 j2Hl, due to the existence of the G3 
flux through Myers' effect |3*7|. Therefore, the junctions of the confining strings 
are dual to the D3-branes, with the electric flux on them embedded in the IR 
geometry. 

The magnetic flux on D3-branes carries the charge of D-strings, which in the 
case of the KS solution have recently been identified with axionic strings 14 7j . 
We do not consider them in this paper. 

The action of a D3-brane in the IR geometry described by the metric (j3J) and 
the R-R flux © is 



S = —T\ 



D3 



J cZVW-det(# a6 + |^F o6 ) +2tt J F 2 AC 2 , (a, b = 0, 1, 2, 3) (6) 



where T str = l/(27ra') and Td3 = l/((27r) 3 a' 2 g stT ) are the tensions of the fun- 
damental strings and D3-branes, respectively. Also, C 2 is the R-R two-form 
potential and G3 = dC 2 . When a confining string consisting of n elementary 
strings is realized as a D3-brane configuration with the topology R 2 x S 2 , the 



4 



n 



integer n is defined as the fundamental string charge carried by the D3-brane. 
The fundamental string current is derived from the action © by differentiating 
it with respect to the NS-NS two-form potential B 2 . (In the derivation of the 
current, F 2 should be regarded as the -B-gauge invariant field strength dA\ + B 2 . 
Once we have obtained the current, we set B 2 = 0.) n is obtained by integrating 
the current as 

- f x dS a D a + f x C 2 = ~ f x dS a D a + J y G 3 , (a =1,2, 3) (7) 

where X is a non-trivial 2-cycle in the D3-brane worldvolume, and Y is a 3-disk 
with boundary BY = X. The flux density D a in (J2J) is defined by 

Da = hW L - ( fl = W) ( 8 ) 

Note that we use Sbi, the Born-Infeld part of the action (JSJ), to define D a . The flux 
density defined in this way is invariant under gauge transformations of the R-R 
2-form potential C 2 . Due to the ambiguity in the choice of Y in S 3 , n is defined 
only mod N. n is a conserved charge and does not change under continuous 
deformations of the brane. 

We define a "flux angle" 6 )f by 

' = < 9 > 

for later use. We fix the mod-7r ambiguity of 9 { by stipulating that < 9 f < ir. 
Because the sum of charges n of the strings meeting at a baryon vertex must be a 
multiple of N, the sum of the corresponding flux angles is a multiple of tt. If there 
are three branches, the sum can be 7r or 2n. These two are essentially the same, 
due to the charge conjugation 6 { — ► it — 6 { . If there are more than three branches 
there are more than two genuinely different cases, as we see below explicitly for 
four-string junctions. 

To clarify the dependence of the action © on the parameters R and N, 
we factor these parameters out of the induced metric and the R-R 3-form flux, 
writing them as 

9ab = R 2 9ab, G 3 = 7^^3, (10) 

where uj 3 is the volume form of S 3 normalized so that its integral over S 3 gives 
$33 = 27r 2 , the volume of the unit 3-sphere. Next, we introduce a rescaled field 
strength F 2 and the corresponding potential A\ as 

* r F2) Al = fL^L Av (11) 



2tt " x 2vr 
The D3-brane action rewritten in terms of these rescaled fields is 



S = - 



r I ~ — r ~ NT. 

Jd i ayj-det( < g ab + F ab )+pjA 1 Au 3 , P= ■ ( 12 ) 
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Here, we have changed the normalization of the action from S to S, which are 
related by S = Td 3 R 4 S. This does not affect the classical equations of motion in 
which we are interested. The dimensionless quantity p is the unique parameter of 
this rescaled action. It is related to the parameter 6 used in Ref. [36J by p = 2/b. 
Also note the relation 

~ = b = ^ ( 13 ) 

where Rb5 = (Na'gstr) 1 ^ 2 is the radius of S 3 near the horizon of N coincident 
flat D5-branes in the flat Minkowski background. Because the MN solution can 
be regarded as the near horizon geometry of D5-branes wrapped around S 2 , 
the radius R is identical to Rr>5, and 6=1. However, for the KS solution, 6 
is modified due to the existence of anti D5-branes, and its numerical value is 
6 = 0.932660368 • • [JTJ EH]- From the action (|12j). it is seen that we can regard 
the dimensionless parameter p as a charge density coupled to the gauge field A±. 

In this paper, we only consider static configurations, and we always omit the 
time dimension in the following. By the term "worldvolume" we refer to a time 
slice of the entire worldvolume. 

Instead of D a defined by (jBJ), we use the rescaled flux density and its Hodge 
dual defined by 

D a = ^ = %^D a , D 2 = l^fdetg ab e abc D a da b Ada c . (a, 6, c = 1, 2, 3) 

(14) 

The relation (J7J) rewritten in terms of the flux angle and the rescaled flux density 
becomes 

27vp6 i = - j D 2 + p j u z . (15) 
This is the integral form of the Gauss's Law constraint, 

dD 2 = P^lbrane, (16) 

where cc>3|brane is the pull-back of U3 to the D3-brane worldvolume. 

The energy of a D3-brane is given by E = T^R 3 E where E is the rescaled 
energy 

E = J rfV^det g ab \J 1 + g ab D a D b . (a, 6 =1,2, 3) (17) 

The shape of a D3-brane and the electric flux density on it should be determined 
so that energy is minimized. 

In the following sections, we use only rescaled quantities, such as E, D a and 
g ab . In the rest of this section, we summarize the relations between these dimen- 
sionless quantities and dimensionful quantities in boundary field theories. From 
the definition of the rescaled quantities S and g ab , we can determine the rela- 
tions between the rescaled energies and tensions and original ones. To determine 
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the corresponding quantities in boundary field theories, we should also take ac- 
count of the warp factor W, which determines the ratio of the energies in the IR 
geometry and those in the boundary field theories. Doing this, we obtain 



(27 



Vertex - WT m R 6 E - ) Jf V ^vertex, 



*™ = W>T m #f=(-^j^Lf, (19) 

where -E V crtex and T are the energy of a baryon vertex and a confining string ten- 
sion, given in the following sections, and E^^lx and T^ FT ^ are the corresponding 
quantities in the boundary field theories. 

The confinement scale Aqcd in the boundary gauge theories can be defined 
by T( FT ) ~ Aq CD . Thus, up to a numerical factor, we have 

W 

Aqcd ~ ^ (20) 

(Here we have used the tension of an elementary confining string, T ~ A r_1 .) 
Then, the energy of a baryon vertex can be rewritten as 

^vertex ~ AqcD^ 3 ^ 2 ^ 2 -^vertex- (21) 

This depends on the string coupling constant, g stT . In the KS case, g str is related 
to the parameter g by (JHJ), and we obtain the following relation, which includes 
only field theory variables: 

42 ~ A QCD iV 3 / 2 ^ vcrtcx . (22) 

By contrast,, there is no parameter corresponding to g str in the MN case. We 
return to this problem in ^ED 

In the following sections, we use only dimensionless quantities, and we omit 
the tildes on rescaled variables for simplicity. 



3 Confining strings and their junctions 

In this section we first briefly review how the tensions of confining strings are 
computed as the energies of D3-branes, following Ref. jHE]- Next, we study the 
kinematics of junctions by considering the balance of tensions on vertices. 

Let (x, y, z) be the three spatial coordinates of the boundary. Combining 
these with the coordinates of the internal space S 3 , we have the set of coordinates 
(x, y, z, 6, <p, ip) with the following metric: 

ds 2 = dx 2 + dy 2 + dz 2 + d9 2 + sin 2 6(dtf + sin 2 0# 2 ). (23) 
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The ranges of the angular coordinates are 



0<9 < 7T, 0<<f)<7i, 0<^<2tt. (24) 

With this parameterization, the volume form of S 3 is given by 

lu 3 = sin 2 9 sin <pd9 A d(p A dtp. (25) 

The dual configuration of an infinitely long confining string is a D3-brane with 
the worldvolume S 2 x R. We begin with the ansatz 

x = a\ y = 0, 2 = 0, 9 = 9\ <p = a 2 , ?p = a 3 . (26) 

The parameter 9 r is a constant representing the angular radius of S 2 C S 3 . By 
rotational symmetry, we can easily determine the flux density D a as a function 
of 9 f and 9 r from (|15|). Its only non- vanishing component is D x , and it is given 
explicitly by 

D x = ~l \9 l -(9 r -sm9 r cos9 r )]. (27) 
b sin 9 r 

Substituting this into (fTTjl . we obtain the tension of a confining string. 



T = ^ = 47rWsin 4 9" + ^(9 l - 9" + sin 9 r cos 9 1 ) 2 . (28) 
dx V 

We should determine the angle 9 r as the point of minimum tension. The condition 
dT/d9 T = gives the relation between # f and 9 r , 

9 1 ' - ( 1 - b 2 ) sin 9 r cos 9 r = 9 { . (29) 

If < b 2 < 1, this relation defines a one to one mapping between < 9 f < it and 
< 9 T < 7T. The minimum value of the tension is 



T = 47rsinfl r Vsin 2 # r + & 2 cos 2 # r . (30) 

For the MN solution, (|2l?jl can be solved immediately, and the tension (jHOJ) 
reduces to the simple form 

e f = 9 T , T = 4vr sin fl f , for b = 1. (31) 

Contrastingly, the angle 9 r and the tension for the KS solution can only be ob- 
tained numerically. The tensions T and angles 9 r for several values of # f are listed 
in Table d These values are used below to compute the energies of the baryon 
vertices. 

In the cases of both the MN and KS backgrounds, the tension depends non- 
linearly on 9 { . This non-linear dependence implies the formation of truly bound 
states and the absence of supersymmetry. Therefore we cannot apply the method 
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Table 1: The angle 9 r and the tension T for several values of 9 
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2tt/12 
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0.785398 
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9.168758 


4tt/12 
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1.047198 
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11.047175 


5tt/12 


= 1.308997 


1.308997 


12.138181 


1.338189 


12.185588 


6tt/12 


= 1.570796 


1.570796 


12.566371 


1.570796 


12.566371 




Figure 2: The balance of three tensions. 

employing Killing spinors, which is useful to determine the baryon configuration 
in M = 4 Yang-Mills theory [T%]. to configurations investigated in this paper. 

Using the tension formula obtained above, we now consider the kinematics 
of three- and four-string junctions. We first study three-string junctions. The 
angles at which the strings meet are fixed by the requirement that of the tensions 
balance at the vertex. If we define the angle 9\ as the angle opposite the force 
vector of the i-th string in the triangle consisting of three force vectors, the angle 
made by the 2-th and j'-th strings is d\ + 0* (see Fig. [2J). If the tensions of the 
three strings are Tj, i = 1, 2, 3, the angles are uniquely determined by 

rp2 I rp2 rp2 

e\ = cos- 1 tiJ-tz — £i o < e\ < 7T, (32) 

i 1 V ) 

and similar equations for 9\ and 9\. These are equivalent to the following: 

T x : T 2 : T 3 = sin 0* : sin 9\ : sin 0* , 9\ + 9\ + 9\ = tt. (33) 

Unlike the relation between 9\ and 0j, 9\ is not a function of only the single 
ang le 9\ with the same index i. Rather, it is a function of the set of angles 
{9{, 6\i 9\}. Because the sum of the three flux angles is ir or 2tt, only two of them 
are independent. We represent 9\ as a function of 9\ and 9\ and denote it by 
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6\ = t(02,OD- The other two angles are similarly given by Q\ = t{6\,6\) and 
0* = r(6[, 91), with the same function r. 

For the MN case, (|33J) reduces to the simple relation sin 9\ = sin 6l Even in 
this case, we cannot determine e\ from only the single angle 6\, because there are 
two solutions, 6\ = 6\ and e\ = n — 6\. These two solutions correspond to the 
two possible values 7r and 2ir of the total flux angle, 0{ ot = J2i=i e\, respectively. 
These two cases can be expressed together by the single equation 

e\ = T{ele f 2 ) = \<K-e{-el\. (34) 

We now proceed to consider four-string junctions. We assume that the four 
external strings are on a plane. We label the external strings in counterclockwise 
order by 1, 2, 3 and 4, and denote the angle made by strings i and i + 1 by 0^+1. 
[see Fig. El (a)] (The labels of the strings are defined mod 4.) The sum of these 



small k, large k 




(a) (b) 

Figure 3: A junction with four external strings. 

four angles is 2-7T, and only three of them are independent. If the flux angles 
{e{, 02, el, el} of the four branches are given, the tensions are determined by the 
formula above, and the balance of these four tensions imposes two conditions on 
the angles 0^+1- These conditions leave one degree of freedom unfixed. We define 
a parameter k parameterizing this unfixed degree of freedom. One useful choice 
is 

« = 013 + 034 = 27T- (023 + 041). (35) 

If this parameter changes, the junction deforms as in Fig. Efb). 

There are three possible topologies of four-string junctions, as shown in Fig. 
@] For junctions with the topologies (H) and (I), which include two three-string 
vertices, the angles 0^+1 are uniquely determined by the flux angles, and the 
corresponding ^-parameters are given by 

k h = 2vr - r{e[, el) - r{e{, e{), Kl = r(e { 2 , el) + r(e{, e[). (36) 

Let us suppose that there is a junction with the topology (H). This can be 
in equilibrium if k = kh- If we try to change the parameter k by changing the 
directions of the strings, the two vertices move in such a way that the movement 
compensates for the variation of k. If we attempt to decrease k, two vertices move 



10 




4 14 14 



(H) (I) (X) 

Figure 4: The three possible topologies of four-string junctions. Two of them 
include two three-string vertices and the other has one four-string vertex. 

away from each other, with the result that k remains unchanged. If we move the 
strings in the opposite direction, the two vertices approach each other. This 
compensating mechanism consisting of the movement of vertices is effective until 
the two vertices coincide. If we continue to move the strings after the vertices 
come to coincide, the topology of the junction tends to change to (I) via (X). 
Until the instant at which the two vertices first coincide, we have k = Kh, and 
the behavior of the junction after this instant depends on the relation between 
kh and Kj. 

If kh > ki, just as the two vertices of (H) meet, the topology of the junction 
changes to (I) and becomes unstable, and then the two vertices separate rapidly, 
eventually realizing the equilibrium condition, k = kj. In this case, the junctions 
of the topology (X) are unstable. 

If kh < Ki, even after the two vertices coincide, the topology does not change 
to (I), and (X) is stable until k reaches Kj. When k reaches ki, the topology of 
the junction finally changes to (I). 

When kh = ki, there is a unique value of k that gives stable configurations, 
and the three topologies can change to each other marginally without changing 

K. 

To express these types of behavior of the junctions, it is convenient to de- 
fine the parameter A = kh — kj representing the "repulsive force" between two 
vertices. If A is positive, two three-string vertices repel each other, and they 
cannot merge to form a stable four-string vertex. If A is negative, two vertices 
attract each other, and they are bound into a four-string vertex, provided that 

Kh < k < Kj. 

Which cases are realized for junctions in the MN and KS backgrounds? From 
()35|) and (j2HJ), A = k h — kj is given by 

A = 2tt - t(9[, 91) - T (6£, 61) - r(9l, 6{) - r(& 9[). (37) 

For the MN solution, the explicit form of A is 

a = 2tt - |tt - e[ - e\\ - |tt - #2 - 03 1 - k - 03 - e{\ -\n-el- e[\. (38) 

We can easily show that A = when J2i=i ®\ — n or J2i=i ®\ — 3vr and that A > 
when £-=i e\ = 2ix. 
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For the KS solution, although it is difficult to determine the signature of A 
analytically, we can numerically check that it is always positive. This implies 
that planar four-string junctions are always unstable in the KS background. 

4 Reduction to lower dimensions 

In the last section, we saw that junctions are described as three-dimensional 
surfaces in the six-dimensional space spanned by the coordinates (x, y, z, 9, 0, -0). 
We only consider planar junctions and set z = 0. If we ignore the extension along 
the (x, ?/)-plane, the worldvolume of a D3-brane representing a fc-string junction 
is S 3 with k punctures. Each puncture is topologically a three-dimensional disc 
and corresponds to each branch string. Let us assume that these punctures are 
located on a large circle of S 3 . In this case, the configuration possesses U(l) 
symmetry, as determined by the symmetry group of this circle. It is convenient 
to choose S 1 given by g^p = as the fixed large circle so that the U(l) action 
constitutes a constant shift of the coordinate ip. General brane configurations 
with this U(l) symmetry are 

x = x(a\a 2 ), y = y(a 1 ,a 2 ), 9 = 9(a\a 2 ), cf> = <f>(a\ a 2 ) , ^ = a 3 . (39) 

Thus, we can represent configurations as two-dimensional surfaces in the four- 
dimensional space with the coordinates (x,y,9,<f>). 

Instead of the angular coordinates (9,(f>), it is convenient to use (u, v,w) sat- 
isfying constraints u 2 + v 2 + w 2 = 1 and w > 0. The relations among these 
coordinates are 

u = sin#cos0, v = cos 9, w = sin 9 sin <fi (40) 

(see Fig. |5J). If we use (x,y,u,v) as a set of independent coordinates, the target 
space is R 2 x D 2 , where R 2 is spanned by x and y, and D 2 is the unit disc sat- 
isfying u 2 + v 2 < 1 in the (u, v )-plane. The variable w = \/l — u 2 — v 2 is treated 
as a function of u and v. Brane configurations are two-dimensional surfaces em- 
bedded in this four-dimensional target space. For example, a confining string 
solution given in is represented as a band in the four- dimensional space which 
is a direct product of a chord v = cos9 r on the unit disc and an infinitely long 
line in the (x, ?/)-plane. 

We should redefine the electric flux density as a field on the two-dimensional 
surface by integrating the flux density D2 over ip. We introduce a one-form D[ 
and its dual vector D' a as 

D[ = ( D 2 , D[ = Jdetg ab e ab D' a da b . (a,b= 1,2) (41) 
J tp 

The Gauss's Law constraint satisfied by D[ is obtained by integrating the con- 
straint (|16|) given by 

dD[ = 2-npdu A dv. (42) 
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Figure 5: The projection of the hemisphere parameterized by (9, <ft) onto the 
unit disc in the (u, w)-plane. The cross-section of a confining string solution is 
represented as a chord of this unit circle in the (u, v )-plane. 



The corresponding integral form is 

2np9 l = - j> D[ + 2np J du A dv, (43) 

where C and A correspond to X and Y in (|15jl. respectively. C is a curve on 
the surface connecting two boundaries. A is a two-dimensional surface whose 
boundary consists of C and a curve in the boundary u 2 + i; 2 = 1 of the four- 
dimensional target space. To obtain the right-hand sides of Eqs. f!42|) and (j43jl . 
we have used the relation 



/ u; 3 = 2n sin 2 sin A d(p = 2irdu A dv. (44) 
The energy (fTTjl integrated over ^ is 

E = J dV/^/HH^. (a, 6 = 1,2) (45) 

Now, the problem that we wish to solve is to find a two-dimensional surface in 
the four- dimensional space (x, y, u, v) and a flux density D[ on it which minimize 
the energy (|45|) under the constraint ()43|) imposed on each external confining 
string. 

When we perform a numerical computation to determine a brane configu- 
ration, we start from an initial configuration and look for a configuration with 
minimum energy by varying the shape of the surface and the flux on it. We 
can choose any configuration as the initial configuration, as long as it satisfies 
appropriate boundary conditions and ()43j) . and as long as it is in the same topo- 
logical class as the final configuration we wish to obtain. One natural choice 
for a confining string is a D3-brane with a shrunk worldvolume, for which the 
second term of ()15j) vanishes. This represents a brane before being expanded by 
Myers' effect. This choice, however, is very singular and not appropriate for the 
numerical computation. Therefore, we instead adopt a worldvolume on which 
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the electric flux density vanishes everywhere on the surface. In this case, the 
constraint (|43|) requires that the area in the (u, t>)-plane enclosed by C and the 
unit circle u 2 + v 2 = 1 be the same as the flux angle 9 { . The most convenient one 
satisfying this condition is the "wedge configuration" defined as a direct product 
of a line in the (x, ?/)-plane and a wedge in the (u, w)-plane, as shown in Fig. El 
(a). 




(a) (b) (c) 

Figure 6: The projections of brane configurations in the (u, t>)-plane. (a) The 
wedge configuration used as an initial configuration, (b), (c) The confining string 
solutions in the MN and KS backgrounds, respectively. 

Starting from a wedge configuration with a central angle 29 { , we can numeri- 
cally reproduce the corresponding confining string solution represented as a chord 
in the (u, u)-plane. Its length is determined by 8 r defined in For the MN so- 
lution, the two angles 8 f and 8 T are equal. This implies that the distance between 
the two endpoints of the wedge on the unit circle is the same as the length of the 
chord [see Fig. ED(b)]. However, this distance becomes slightly greater for the KS 
solution, [see Fig. ED(c)]. 

Wedge configurations are also suitable to construct initial configurations for 
junctions, because they can be pasted just like interaction vertices in Witten's 
open string field theory as shown in Fig. [7| (a). By varying this initial 




(a) (b) (c) 

Figure 7: (a) An initial configuration for a three string junction is composed of 
three wedge configurations. After minimizing the energy by varying the configu- 
ration, we obtain the junction solutions for (b) the MN background and (c) the 
KS background, respectively. 

configuration and minimizing the energy, we obtain a junction configuration. 

In the case of the MN background, the distance between the endpoints of each 
curve representing each branch string does not change, and the junction solution 
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is expected to be a triangle inscribed in the unit circle in the (u, t>)-plane [see Fig. 
13(b)]. Indeed, the coincidence of the endpoints of the three chords is guaranteed 
by the relations (|53|) . given below. It is important that this triangle is similar to 
the tension triangle in Fig. El 

For the KS solution, the chords representing asymptotic string solutions are 
too long to make an inscribed triangle [see Fig. [7| (c)]. 

5 Numerical analysis 

In this section, we report our results of our numerical study of brane configu- 
rations. We start by reconstructing the confining string solutions to assess the 
accuracy of the results of the numerical method by comparing them to the ana- 
lytic results given in JJJ 

We first prepare a wedge configuration with a central angle 2# f and a length 
L as an initial configuration [see Fig. |S] (a)]. We compute the energy of the 
surface by means of the triangulation method. The surface is divided into a 



Figure 8: (a) An initial wedge configuration and (b) a numerically generated 
confining string solution. The segments in the worldvolume in (b) represent the 
flux density. 

mesh of 2n^ csh small triangular regions. The shape of the surface is represented 
by the positions of the sites, and the flux density is treated as link variables. 
We seek a configuration that minimizes the energy by varying these variables 
while enforcing the Gauss's Law constraint. The detailed algorithm is given in 
the appendix. A result for {9 l = n/3, L = 4, b = 1) is given in Table El By 
comparing this result to the corresponding tension in Table ^ we can estimate 
the accuracy of the outputs of the numerical analysis. We find a relative error 
of order ~ 10 -4 for the finest meshes, with n mes h ~ 80. Plotting the error as a 
function of the typical size of triangles, a = n~ esh , we find that the error is almost 
proportional to a (see Fig. EJ). The value for n mes h = oo (a = 0) given in Table 
El is obtained by extrapolating the data for finite n mcs h with a polynomial of a. 
This extrapolation gives a value very close to the exact one, with a relative error 




(a) 



(b) 
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Table 2: A numerical result for the tension of a confining string with Q l = 71/ 3 in 
the MN solution. The value for n mcs h — > 00 is obtained by extrapolating the data 
for n mcsh = 10, 20, 40 and 80 with the fitting function c + Cin~g Sh + c 2 n~g sh + 

-6 
c 3^mesh- 



^mesh 


10 


20 


40 


80 


00 


exact 


T = E/L 


10.7858 


10.85978 


10.87722 


10.88144 


10.88283 


10.88280 


error 


-0.0970 


-0.02302 


-0.00558 


-0.00136 


+0.00003 






Figure 9: A fitting of numerically generated tensions for several values n mes h. 



< 10~ 5 . Because of this small error, we use the same extrapolation to obtain the 
energies of the baryon vertices. 

Let Ej unc (L) be the energy of a three string junction with the branch length 
L. Then, the contribution from a baryon vertex is defined by 

E veitex = lim (£ junc (^ f , L) - L J2 T(6{)) . (46) 

I 

However, we cannot in practice numerically compute E }vcac in the limit L — > 00. 
Therefore, instead, we use a sufficiently long fixed length beyond which the effect 
of the vertices becomes negligible. In the following way, we find evidence that 
L = 4 is sufficiently large to to guarantee a relative accuracy of 10 -5 . First, we 
compute the energies of two strings stretched between x = and x = L with 
lengths L = 4 and L = 8. The worldvolume of the short and the long strings 
are composed of 2n^ esh and 4n^ esh triangles, respectively. Hence the areas of the 
triangles are almost the same. The difference between the present computation 
and that described above to estimate the accuracy is that here on one boundary 
at x = we impose a fixed boundary condition that prohibits the shape of the 
boundary from changing from a wedge to a straight chord. Thus, the resulting 
configurations have a wedge-shaped boundary at x — 0, and as x becomes large, 
the shape approaches a string solution (see Fig. ITU]) . If the difference in energy 
between the two branes of lengths L = 4 and 8 is sufficiently close to the energy of 
a confining string of length 4 = 8 — 4, we judge the configuration to be sufficiently 
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Figure 10: The brane configuration with a boundary on which the fixed boundary 
condition is imposed. 

close to the string solution around x ~ 4, and we can use L = 4 to compute the 
energies of the baryon vertices. The result of this analysis is given in Table 

Table 3: The numerically computed energies of the confining string configurations 
with a fixed boundary at one end. T(n mes h) is the numerically computed tension 
given in Table El 



^mesh 


10 


20 


40 


80 


El=8 


87.2002 


87.75829 


87.88710 


87.91770 


El=4 


44.0562 


44.31901 


44.37818 


44.39194 


(E L=8 -E L=i )/(8-A) 


10.7860 


10.85982 


10.87723 


10.88144 


Difference from T(n mes h) 


+0.0002 


+0.00004 


+0.00001 


+0.00000 



El There it is seen that indeed the configuration is sufficiently close to a string 
solution before x reaches L — 4. 

Based on the above preliminary results, we can finally compute the energies of 
the baryon vertices. The initial configurations are combinations of three wedges 
pasted like Witten's 3-string vertices [see Fig. [TT] (a)]. The directions of three 
branches of initial configurations in the (x, y)-plane are set as in Fig. EJ so that 
the tensions are balanced. Now, instead of taking the L — > oo limit, we simply 
use L = 4, and each branch consists of 2n^ esh triangles. 

We give the results for the KS background in Table |U Because vertices with 
more than three branches are unstable, we only give the results for three string 
vertices. Contrary to the intuitive expectation, the signatures of the baryon 
vertex energies are negative. Their absolute values depend on the ratio of the 
flux angles, and it seems that the energy becomes closer to zero as the flux angles 
become increasingly asymmetric. This behavior was expected because in the 
limit that one of the flux angles vanishes, the junction becomes a confining string 
configuration, and the energy of the vertex goes to zero. 
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(a) (b) 

Figure 11: An initial configuration composed of three wedge configurations and 
the corresponding final junction configuration obtained numerically. 

Table 4: The energies of baryon vertices in the KS solution. Ej nnc is the energy 
of an entire junction of branch length L = 4. Eq represents the contribution 
of branches, computed as E = Z)f =1 LT{9\). The energy of a vertex is defined 
by -Evcrtcx = -Ejunc — Eq. The shaded numbers are the final results obtained by 
extrapolation. 
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123.8908 
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-1.418 


-0.5559 


-0.3488 


-0.2981 


-0.2813 



The results for the vertex energies in the MN background are more interesting. 
As shown in Table El the ratios of the vertex contributions to the total energies 
are smaller than 10 -5 . This strongly suggests that the energies of baryon vertices 
exactly vanish for any values of the flux angles. In other words, the energy of 
a junction is obtained by summing up the contribution of each branch given as 
the product of its tension and its length in the (x, y)-plane. In Table El we give 
another result for a stable planar four-string vertex. It too seems consistent with 
a vanishing vertex energy. 

Before ending this section, we confirm that junction configurations in the MN 
background have the asymptotic forms suggested in which are represented as 
triangles in the (w, w)-plane. The projection of a numerically generated junction 
in the (u, t?)-plane is shown in Fig. El and it indeed seems to be a triangle, as 
expected. 
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Table 5: The energies of baryon vertices in the MN solution. Ej uncj Eq, and 
-^vertex are defined in the same way as those in Table EJ The shaded numbers are 
the final results obtained by extrapolation. 
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/ 2n bn bn \ rp 
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Table 6: The energy of a planar four-string vertex with 9\ = 7r/4. Ej unc , Eq, and 
EVertex are defined in the same way as those in Tabled and The shaded number 
is the final result obtained by extrapolation. 
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6 Analytic proof 

As we saw in the last section, the numerical results for junctions in the MN back- 
ground strongly suggest that baryon vertices do not contribute to the energies 
of junctions. This is also the case for supersymmetric junctions consisting of 
(p, g)-strings. [HJ E01 EH E3 ESI EH EH] This fact for (p, g)-junctions is proven by 
dualizing junctions to membrane configurations in M-theorv|54j. From this view- 
point, junctions are represented as two-dimensional smooth surfaces embedded in 
a four- dimensional space R 2 x T 2 . Let (X, Y) and (U, V) be orthogonal coordi- 
nates in R 2 and T 2 , respectively. We combine them into two complex coordinates, 
Z\ = X + iU and Z2 = Y + iV. Then, a surface in this space is represented by 
a function Z2 = F(Z\, Z\). It is known that if we require a brane configuration 
to be supersymmetric, the function F must be holomorphicjHU EHj- This fact is 
very important in proving the disappearance of the energies of vertices [5*4*]. 

Brane configurations in the MN background have a similar property. As 
we have seen, they are described as two-dimensional smooth open surfaces in 
the four-dimensional space R 2 x D 2 with coordinates (x, y, u,v). Let us define 
complex coordinates Z\ — x + iu and z 2 = y + iv in analogy to the (p, g)-junction 
case. A surface can be represented by an equation z 2 = f(z\,Zx). Far from the 
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Figure 12: A three-string junction configuration in the MN background projected 
on the (u, v )-plane. 

vertex, each branch of the surface asymptotically approaches a confining string 
solution, which is represented by 

y = diX + bi, u = CiV + d iy (47) 

where i is the label of the branch. When we rewrite these relations in terms 
of complex variables z\ and Z2, the function / is holomorphic if the coefficients 
Oj and Cj are the same. We can always realize a, = c, for one i by applying 
the appropriate rotation of the coordinates u and v. Furthermore, due to the 
similarity of the triangles displayed in Fig. El and Fig. [7[ the relation Oj = 
Cj can be realized for all % simultaneously. Therefore, the function / becomes 
holomorphic in the asymptotic part of all the branches. We call this property 
"asymptotic holomorphy". In terms of real variables, the holomorphy of the 
function / is represented by the Cauchy-Riemann relations: 




In fact, the second equation in (|4*H|) is vacuous as a condition for the asymp- 
totic shapes of surfaces. Because we are considering brane configurations which 
asymptotically approach confining string solutions, the asymptotic forms are al- 
ways represented by the two factorized equations in (|47|). Hence, here, the second 
equation in PK|) reduces to the trivial relation = 0. 

This asymptotic holomorphy can be used to check the stability of four-string 
vertices. Under the assumption of the coincidence of the endpoints of chords, 
the relation { = 9 t for the MN background implies that k chords representing 
branches of a /c-string junction form a fc-gon inscribed in the unit circle in the 
(u, w)-plane. For four-string junctions, there are three cases distinguished by their 
values of 9{ ot = Ya=\ If #tot = 71 or 3vr, the squares defined in the (u, v) -plane 
in this way are similar to the squares composed of the four tension vectors of 
four external strings, just as in the case of three-string junctions. However, if 
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6{ ot = 2ir, the inscribed squares are twisted and cannot be similar to the tension 
squares [see Fig. (&)]• in this case, we obtain squares similar to the tension 




v 



u 



(a) 



(b) 



Figure 13: (a) A twisted square for 6{ ot = 2tx. (b) Two ways of unfolding. 

squares by unfolding the twisted squares. The two ways of unfolding correspond 
to the two different topologies (H) and (I), which cannot be deformed smoothly 
into each other [see Fig. H3l(b)]. From these facts, we conclude that a four-string 
vertex is stable only if it satisfies the asymptotic holomorphy condition, which is 
equivalent to the similarity of the tension squares and the inscribed squares in 
the (u, v )-plane. 

The asymptotic holomorphy relations are not sufficient to compute the en- 
ergies of baryon vertices. For this purpose, we need to find some relations that 
hold everywhere on branes. One possibility is to extend the asymptotic holomor- 
phy to the entire worldvolume as they are. Unfortunately, however, a numerical 
investigation of generated junction solutions demonstrates that this cannot be 
done. 

In the case of super symmetric configurations, there is a technique to deter- 
mine relations holding everywhere on the surface by using Killing spinors. In- 
deed, various supersymmetric embeddings of D5-branes in the MN background 
are obtained in Ref. ^2] with this technique. Unfortunately, this is not appli- 
cable to junction configurations, because they break all the supersymmetries, as 
mentioned in JI1 

The absence of systematic methods compels us to resort to guesswork. Nu- 
merically generated junction configurations give important information about 
how to approach to this problem. Let us consider Fig. The radial pattern of 
the flux in it implies that (D U ,D V ) oc (u, v). The proportionality factor can be 
determined by the Gauss's Law constraint (|42|). and we obtain the ansatz 



where u a = (u, v) and we use the static coordinates a a = u a . Then, by substitut- 
ing this ansatz into the energy (|45|). we obtain 



D[ = 2ir(udv - vdu), D 



Vdet g ab 



u 



a 



(49) 




(50) 
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with the effective induced metric g a b defined by 

g ab = -Sab + wX l a X l b , (51) 
w 

where X 1 = (x,y) and X l a = dX l /du a . The expression (J50|) implies that the 
two-dimensional surfaces can be treated as Nambu-Goto-type branes in the back- 
ground with the effective metric 

— 2 1 

ds = — (du 2 + dv 2 ) + w(dx 2 + dy 2 ). (52) 
w 

Thus, it is quite natural to replace du, dv, dx, and dy in the Cauchy-Riemann 
relations (jlHj) by w~ x l 2 du, w~ x l 2 dv, w l / 2 dx, and w l l 2 dy, respectively. As a result, 
we obtain the following "modified Cauchy-Riemann relations", which take the 
place of the holomorphy for (p, g)-string junctions: 

(53) 

Note that these relations do not contradict the asymptotic holomorphy relations 
(EHJ) but, rather, guarantee them. Because the relation between (x, u) and (y, v) 
is factorized as ()47|) in asymptotic regions, the second relation in (p)3*j) becomes 
trivial as the second relation in (|48j). 

We can show that the equations of motion are automatically satisfied if the 
shape of a brane satisfies the modified Cauchy-Riemann relations (|53|) and the 
flux density on it is given by the ansatz (|49|) . The equation of motion for X % 
obtained from the energy (joHj) is 






d a ^det g ab wg ab XI j =0. (54) 

In order to obtain the equation of motion for the gauge field, we should use the 
action (@EJ), to which the ansatz ()49)1 has not been applied. With variations in the 
form 5D' a = (l/^/g)e ab db4>, which do not violate the constraint (|4"2*|) . we obtain 
the Maxwell's equation for the electric field strength E a , 

e ab d a E b = 0, (55) 

where the relation between E a and D a is 

E a = , 9abDa (56) 
'(2irw) 2 + g a bD a D b 



This can be rewritten using the ansatz ()49j) as 



p 9abU b g a bU b . , 

E a = = — = . 57 

Vdet g ab uy det 9ab 
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To show that equations of motion (|54|) and (|55|) hold, we first rewrite the 
relations (|53|) in terms of the functions x(u, v) and y(u, v). Then, using the chain 
rule for partial derivatives, we obtain 



•Eu Vv 0, x u y v Uu^v ni (5 



w 2 



or equivalently, 



trX: = 0, detX: = ^. (59) 

(Here and hereafter, when subscripts are used to represent derivatives, the inde- 
pendent variables are always u and v.) Using these relations, the effective metric 
floTj) becomes 



g ab = w{XlXl + 5 ab det X* a ) = w(x v - y u ) yu yv , (60) 

Y Xu X v J 

and its determinant is 

detg ab = (x v -y u ) 2 . (61) 

With the help of (J60|) and (J61j) . we can easily show that the equations of motion 
(1541) and ()55|) are indeed satisfied. 

Although we have shown that the modified Cauchy-Riemann relations (|53|) 
and the ansatz (|49|) guarantee that brane configurations satisfy the equations 
of motion (J54)) and (|55|). the converse is not true, as there exist solutions of 
equations of motion which do not satisfy the modified Cauchy-Riemann relations 
and the D-field ansatz. An example of such a configuration is the deformed 
baryon configuration given in Ref. j^Sl- Such configurations, however, do not 
approach confining string solutions asymptotically, and their existence does not 
invalidate the arguments given in this section. 

Now, we are ready to prove the vanishing of the energies of baryon vertices 
analytically. From (jSOjl and (|61|h the energy Ej unc (L) defined above (pj)) is given 
by 

Ej U nc(L) = 2ir / dudv\x v — y u \ = 2ir <t> (xdu + ydv) , (62) 

J J JdJ 

where J is the region in a junction worldvolume determined by the relation x 2 + 
y 2 < L 2 . Before applying Stokes' theorem to (|81Z|). the absolute value should be 
moved out of the integral. This is possible, because as shown by ([58)1 . x v — y u 
never vanishes on the worldvolume. The boundary dJ consists of two kinds 
of boundaries, which we call A-boundaries (the solid lines in Fig. IT4"|) and B- 
boundaries (the dashed lines in Fig. IT4"|) . The A-boundaries are the edges of 
the surface on the boundary u 2 + v 2 = 1 of the target space. The B-boundaries 
arise due to the cut off x 2 + y 2 = L. The A-boundaries are mapped to apices 
of the triangle when they are projected on the (u, v )-plane, and the coordinates 
u and v are constants on them. Therefore, the A-boundaries do not contribute 
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Figure 14: (a) J and its boundary, (b) The projection on the (x, y)-plane. (c) 
The projection on the (u, -plane. 

to the contour integral in (J62j) . Thus, we need only consider the integral on the 
B-boundaries. If L is sufficiently large, each B-boundary is mapped to a point by 
the projection in the (x, y)-plane. Let (xi,yi) be the point corresponding to the 
z-th B-boundary, Bi. These are on the circle x 2 + y 2 = L 2 , and we can write 

(x ij y i )=L(a i ,/3 i ), (a 2 + [3 2 =l) (63) 

where (ctj, is a unit vector representing the direction of the i-th branch in the 
(x, y)-plane. If we project the B-boundaries to the (u, u)-plane, they are mapped 
to sides of the triangle when L is sufficiently large. Due to the similarity be- 
tween this triangle and the tension triangle, which is guaranteed by the modified 
Cauchy-Riemann relations (JHHJ), the integration measure [du^dvi) on can be 
represented as 

(du h dvi) = ±(ai, Pi)dli, (64) 

where the signature depends on the orientation of the contour, and dli is an 
infinitesimal length on the i-th side of the triangle in the (u, v )-plane. Combining 
(J53J) and (|64[). we obtain the following expression for the energy of a junction: 

E imic = 2ix Y, L I dk = 2nJ2Lx2 S m6l = J2 LT{9{). (65) 

i jBi i i 

Here, we have used the fact that the integral J B . dU gives the length of the z-th 
side of the triangle 2 sin 9\ and the tension formula ()31|) for the MN background. 
The result (J65|) implies that the energy of a junction is given by the sum of the 
energies of the branches and that there is no vertex contribution. 

Now we can answer the question raised at the end of Although the expres- 
sion ()21|) seemingly depends on ^ str , which has no counterpart in the boundary 
field theory of the MN solution, -E vertC x on the right-hand side is actually zero, 
and the energies of the baryon vertices in fact are independent of g stT . 
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7 Conclusions 



In this paper, we have investigated string junctions in TV = 1 supersymmetric 
gauge theories in the context of the gauge/gravity correspondence. We have used 
the Maldacena-Nunez and the Klebanov-Strassler solutions as gravity duals of 
the confining gauge theories. 

In JJ] we studied the balance of the tensions for three-string junctions and 
planar four-string junctions. We found that four-string vertices are stable only 
in the case that 6\ ot = n or 3tt in the MN solution. Planar four-string vertices 
in the KS solution are always unstable. It may be interesting to generalize this 
consideration to non-planar junctions. 

In £0 we reported the results of numerical computations of the energies of 
baryon vertices, and we found that in the KS case the energies are negative, while 
in the MN case they almost vanish. The vanishing of the energies strongly sug- 
gests that the brane configurations in the MN background possess some analytic 
structure which guarantees this vanishing. Indeed, we discovered relations sim- 
ilar to the Cauchy-Riemann relations for holomorphic surfaces. With the help 
of these relations, we analytically proved the disappearance of the energies of 
baryon vertices. 

We should emphasize that the analysis given in this paper is classical. As 
mentioned in JJJ the existence of confining strings breaks all the supersymmetries. 
Thus, there is no mechanism to control quantum corrections, and our results are 
justified only in the limit of large iV and large Ng stT , in which the background 
curvature is small compared to the string scale and the Plank scale. 

In order to investigate realistic baryons in non-supersymmetric QCD, we 
have to use different gravity backgrounds. For example, we can use an AdS 
Schwarzchild black holejU EZ]- Because the IR geometry of this solution has 
structure similar to the KS and MN solutions, we can study junction configura- 
tions in it in a manner similar to that used in this paper. 

We have treated only brane configurations embedded in the IR geometry in 
this paper. They can be used for highly-excited baryons, in which the endpoints 
of strings are separated from each other. To analyze ground-state baryons, we 
should consider different kinds of brane configurations that consist of a baryon 
vertex and only one external string going up along the r direction to flavor branes. 
The external string in this case represents coincident quarks. It is important to 
determine whether such branes are stable, and if so, to investigate their energies 
and excitations. 

There are many interesting problems associated with brane constructions of 
hadrons in addition to those mentioned above. We hope to return to these issues 
in the near future. 
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A 

A.l Treatment of flux on triangulated surfaces 

The purpose of this section is to explain how to descretize branes with flux flowing 
on them and how to vary variables while maintaining the Gauss's Law constraint. 

To compute the energy of a D3-brane, which is represented in (|43|l as an 
integral over a two-dimensional surface, we first triangulate the surface. We label 
sites by i,j, k, . . . and each oriented link by two labels representing its two ends. 
The functions A 7 (cr a ) describing the shape of the surface are replaced by the 
variables X\ = (xi, yi, Ui, Vi, Wj) on the sites, with the constraints uf + vf + wf = 1 
and Wi > 0. To represent the flux density, we use link variables 

(jUj = I' D[, (66) 

J i 

where the integral is carried out along the link ij. This represents the amount 
of flux flowing across the link ij. We define the orientation in such a way that 
if the arrow from site i to j is upward, 0jj represents the flux passing the link ij 
from left to right. By definition we have 

(pij = ~4>ji- (67) 
The area of a triangle ijk is denoted by (see Fig. I15|) . This real 




Figure 15: The triangle ijk. X- denotes the coordinate of the vertex i. <pij 
represents the flux flowing across the link ij. The area of the triangle ijk is 
denoted by s^fc. 
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positive number. We also define s'^ k as the area of the triangle ijk projected on 
the (u, t>)-plane. This take either a positive or negative value, depending on the 
orientation of the triangle, and s' i j k are easily represented as functions of the 
variables X\, Xj and X k . 

Given the variables X\ and 4>ij, the energy of the D3-brane is obtained as the 
discretized version of the energy (^J) . 

E= E s ijk l(2irw ijk y + £(£>g fc ) 2 , (68) 

triangles y / 

where Wij k is the w-component of the center of mass X\- k = (X- + X 1 - +X^)/3 of 
the triangle ijk, and D't k is the push-forward of the electric flux density to the 
five-dimensional space (x, y, u, v, w), which is given by 

D 'm = - x k) + M*j - x m) + M*i ~ x ik)}- (69) 

£S ijk 

A discretized version of the Gauss's Law constraint is 

4>ij + 4>jk + 4>ki = 2nps' ijk . (70) 

In order to find a configuration that minimizes the energy, we should vary the 
variables X\ and 0^- in such manner that does not violate this Gauss's Law 
constraint. 

There are two kinds of variations. Variations of <pij with fixed X\ are generated 
by the following variation for each site i: 

<fiji -> (pji + c, j e Adj(i). (71) 

Here, Adj(z) represents the set of all the sites adjoining the site i. This variation 
changes the rotation of the flux density around the site i [see Fig. (a)]- The 
other kind of variations are those that change X\. Even when we vary the 
variables X-, we should take account of the Gauss's Law constraint (|70|). because 
variations of X\ change charges in triangles. If the position of a site i moves from 
X\ to X\, = X\ + SXf, the projected area s'^ k of the triangle ijk is changed by 
the amount s'j U , — s kii , [see Fig. HEl(b)]. To maintain the Gauss's Law constraint 
(|7()|) . we must change the flux variables simultaneously according to the relation 

<pki> = <Pki - 2nps' kii ,, k e Adj(i). (72) 

Any continuous deformation can be generated by the two kinds of variations (|7T1) 
and 
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Figure 16: (a) A variation changing the rotation of the electric flux density 
around i while maintaining the Gauss's Law constraint, (b) If a site i moves to 
i', the projected areas of the triangles that possesses this site as a corner change. 
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